Unbinding of giant vortices in states of competing order 
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We consider a two-dimensional system with two order parameters, one with 0(2) symmetry and 
one with O(M), near a point in parameter space where they couple to become a single 0(2 + M) 
order. While the 0(2) sector supports vortex excitations, these vortices must somehow disappear 
as the high symmetry point is approached. We develop a variational argument which shows that 
the size of the vortex cores diverges as 1/vA and the Berezinskii-Kosterlitz-Thouless transition 
temperature of the 0(2) order vanishes as l/ln(l/A), where A denotes the distance from the high- 
symmetry point. Our physical picture is confirmed by a renormalization group analysis which gives 
further logarithmic corrections, and demonstrates full symmetry restoration within the cores. 
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Competition between different states of order is a hall- 
mark for a large class of strongly correlated systems. 
Some (quasi-)two dimensional examples are competing 
magnetism and superconductivity in cuprates [T], or- 
ganics [2] and the recently discussed thin-film heavy- 
electron systems [3] ; competing superfluid and crystalline 
order that may occur in bosonic systems on triangular 
lattices [4] or cold atomic gases [6]; or charge den- 
sity wave order that competes with superconductivity in 
Sri4- K Ca a; Cu2404i [7]. In all these cases, at least one of 
the two competing order parameters has XY, i.e. 0(2) 
symmetry, while the other may in general be 0{M). The 
cases M = 1, 2, and 3 correspond to the other order 
parameter being of Ising, XY, or Heisenberg type, re- 
spectively. The M = 1 case describes easy-plane mag- 
netism [SI [S] ; the case M = 2 relates to supersolid phases 
in cold-atom systems [5] and to competing density-wave 
and superconducting order in layered materials [IOHl2] : 
models with M = 3 have been considered in the con- 
text of high-temperature superconductivity [13] . Both 
order parameters interact and the symmetry of the cou- 
pled problem is 0(M) x 0(2). However, at a certain 
fine-tuned point in phase space one may expect the sym- 
metry to be enhanced, from 0(M) x 0(2) to 0(N) with 
N = M + 2. This is not the most general scenario for 
competition between two order parameters, but it has 
been conjectured to occur in many different microscopic 
models, including all of the cases mentioned above [5rH3]. 

This symmetry enhancement acquires a particularly 
interesting aspect in layered or two-dimensional systems, 
where long range order is absent for continuous sym- 
metries due to the Hohenberg-Mermin- Wagner theorem 
P3]. However, the 0(2) sector supports non-trivial topo- 
logical configurations, i.e. vortices. The unbinding of 
vortex-antivortex pairs converts an algebraically ordered 
superfluid or crystal to a disordered normal fluid. This 
Berezinskii-Kosterlitz-Thouless (BKT) transition [T5HT7] 



occurs at a non-zero temperature Tbkt- 

Suppose the fine tuning to a high-symmetry point in 
the phase diagram is achieved by varying a dimension- 
less parameter A > towards A = which corresponds 
to the 0(N) symmetry point. In each realization of 
this model, the experimental handle corresponding to 
our parameter A is different - for example in the con- 
text of cuprates/organics it would correspond to dop- 
ing/pressure |13j while for cold dipolar bosons it may be 
controlled via the angle of a polarizing field [5. . How- 
ever, no matter which particular microscopic realization 
of this model is chosen, if N > 2 then Tbkt must vanish 
for A — > 0. Indeed, combining spin-wave based renor- 
malization group calculations with crossover arguments 
one can estimate that Tbkt vanishes as l/log(l/A) as 
A vanishes - see Fig. [T] This was first derived for the 
case M = 1 in [8]; the present work extends this result 
to generic values of M. There are, however, a number 
of nontrivial aspects that emerge from this picture. The 
BKT transition must vanish because of the dominance of 
spin wave excitations of the high symmetry model. On 
the other hand, spin waves do not usually interfere with 
vortices: in the 0(2) model, spin waves do not lead to 
a renormalization of the stiffness. So how do vortices 
become spin waves? 

In this paper we investigate the fate of vortices of 
XY-order parameters, and of the BKT transition that 
they mediate, as the high-symmetry point is approached. 
Combining variational arguments and a renormalization 
group (RG) analysis, we study the crossover and transi- 
tion temperatures and show that for small A, the size of 
the vortex core diverges as 
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where a is the core size of a single 0(2) order parameter 
(i.e. for A J> 1). This is a consequence of the emergence of 
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competing order and O(N) spin waves inside the vortex 
core. Thus, the enhanced symmetry becomes visible not 
only at high temperatures, where A may be neglected, 
but also at low temperatures, via the size of the vortex 
core near (and below) the unbinding transition. 

In accordance with the usual ideas of universality, we 
consider the long wavelength action with the appropriate 
symmetry, which is a perturbed non-linear sigma model: 



S = 
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Here n = (s, m) is an iV-component vector subject to the 
unit-length constraint n 2 = 1. The vector s has two 
components while m has the remaining M components; 
these two vectors corresponding to the two competing 
order parameters of the original theory. The matrix D is 
given by 
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The model at A = has full O(N) symmetry, whereas 
for A 7^ this symmetry is broken to O(M) x 0(2) by 
giving a mass to the m sector of the theory. The ac- 
tion contains a reference length scale, a (of the order of 
the crystal lattice spacing), and a reference energy scale, 
J. The exact meaning of J depends on the microscopic 
model from which ^ is derived; but typically J corre- 
sponds to a bandwidth of the unperturbed microscopic 
model, possibly reduced by quantum fluctuations or geo- 
metric frustration. We assume we are in a regime where 
the ground state is fully ordered so that any possible 
quantum dynamics beyond ^ may be safely neglected. 

Let us start by considering the situation far from the 
O(N) point, i.e. for sufficiently large values of A. In 
this case the vector n is effectively constrained to lie in 
the easy plane of the A-term, and thus becomes a two- 
component one. Hence the standard BKT argument [151 
fit)] applies here, and the relevant excitations to consider 
are vortices, which are described by a configuration 



(71,1,71,2) = (cos 9, sin t 
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where (r, 9) are plane polar co-ordinates with r = at 
the center of the vortex. The action for a single such 
vortex is 



1 \ a 



S c ,, 



(5) 



where L is the linear size of the system, a is the lattice 
spacing, S^orc {T) ~ it/2T is the action of the (lattice- 
scale) core of the vortex, and we work henceforth in units 
where J = 1. In ^ we were obliged to cut off the di- 
vergent energy near the center of the vortex at the lat- 
tice scale a, because the model contained no other length 
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FIG. 1. [Color online] Schematic phase diagram of the model 
(2} for A > near A = 0. The thick blue line indicates Tbkt 
as a function of A. Below this line the model shows power- 
law correlations; above the line it is XY-disordered. The red 
line shows the crossover temperature between a disordered 
state of 0(2) character (proliferation of vortices) and one of 
O(JV) character (spin waves). The crossover line occurs where 
the size of the vortex core becomes of the same order as the 
average distance between vortices. The arrows show the RG 
flow ( 13 1 for the specific case M = 2 (see text). The shadowed 



must be complemented by vortex considerations. 

scale. As A is reduced, however, a new mechanism of 
removing this divergence becomes available: the vector 
n can simply be rotated away from the easy plane [20] , 
This costs an energy proportional to A, but is worth do- 
ing near the center of the vortex where the vortex action 
would otherwise be very high. Let us, then, consider 
a configuration described by a new length-scale £: for 
r > £, it is identical to the above-described vortex, while 
for r < £, n has components perpendicular to the easy 
plane. The action of such a configuration may be esti- 
mated as 



core 
7T A 

T2^ 



S, 



vortex 
7T 
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minimizing this with respect to £ determines the opti- 
mum core size, 



A 



(7) 



This has the interesting consequence that while the core 
size increases and the total vortex action decreases with 
decreasing A, the core action remains the same: Score = 
7r/2T, as in the original BKT case. 
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The finite-temperature transition from the quasi-long- 
range-ordered to the disordered state occurs via a pro- 
liferation of unbound vortices and anti- vortices. To de- 
scribe this one considers the renormalization group flow 
equations [To] 



d£ 



dy(i) 
d£ 



(2-nT-\l))y(l), (8) 



where T(£) and y{£) are the renormalized temperature 
and vortex fugacity at length scale £ = ae e . These flow 
equations have a separatrix along the line ir/2T = 1 + 
27ry, meaning that the transition temperature satisfies 
the following equation: 



Ibkt = 



2 + 47ry (T BK t) ' 



(9) 



In the regular single scale BKT transition, the bare fu- 
gacity is given by y(T,£ = 0) = e -S ' coro( ' T ) |1B], which 
leads to Tbkt ~ k/2 as the fugacity gives only a small 
correction to this value. As we have seen above, S'corc 
does not depend on A, so we must ask the question: how 
is the BKT transition modified in the presence of another 
length scale £o? This question is much more general than 
the model in this paper, and has recently arisen in dif- 
ferent situations [181119]. 

The RG flow ([§) should now start at the length scale 
£o = o,e e ° and not at the microscopic scale a; so for the 
initial values we need to know the renormalized temper- 
ature and fugacity at this length scale. Studying the 
problem just from the point of view of vortices (i.e. ig- 
noring spin waves), the large core derived above is inert, 
and thus T(£ ) — T(£ = 0). However even for the inert 
core, the fugacity has a naive scaling dimension and thus 
flows according to 



dy/d£ = 2y; 



(10) 



integrating this equation gives y(£o) = {£,o/a) 2 y{£ = 0). 
This enhancement of the fugacity can be understood 
physically by realizing that while the vortices live at a 
length scale £o> the entropy comes from enumerating the 
possible positions for the center of the vortex which in- 
volves the lattice scale a. The equation for Tbkt now 
becomes 



T, 



BKT 
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Substituting in the form |7]) for the optimum core size 
and defining a new variable x = 7t/Tbkt we obtain A = 
Ane~ x / 2 / (x-2). We are interested in the solution of 
this equation as A — > 0, in which case e~ x l 2 / (x — 2) 
must also tend to zero, i.e. x — > oo. Keeping only leading 
order terms, we see that x ~ In (1/A), which gives 



?bkt 
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Thus we see that this simple argument gives a BKT tran- 
sition temperature that vanishes as A — ¥ 0, as expected 
on symmetry grounds. Furthermore, it defines the length 
scale £o : below this length scale the physics of the sys- 
tem becomes sensitive to the proximity to an enhanced- 
symmctry point; above it, the physics is essentially that 
of the large-A, 0(2) system. 

The analysis above is from the point of view of vortices; 
it doesn't include spin-waves. In particular, we should 
take into account that the full O(N) dynamics may still 
be intact inside the vortex core. To analyze this issue we 
turn to the renormalization group treatment of spin- wave 
excitations. The renormalization flow equations of O(N) 
nonlinear sigma models with symmetry broken by giving 
M of the N components a mass have been studied in [2JJ ; 
they also follow from a generalization of the argument for 
the case M = 1 given in [5]. For generic N and M we 
obtain: 



dT(£) _ T{lf 
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where T(£) and A(^) are the renormalized temperature 
and anisotropy at length scale £ = ae l as defined before. 
For the present case N = M + 2, the first flow equation 
simplifies to 



dT{£) T{£) 2 N-2 



dt 



2tt 1 + A(£) 



pa| ) 



As long as A(£) is small, the flow of the temperature (i.e. 
of the inverse stiffness) is that of the usual O(N) model, 
while dT jp — > at large A(£) as expected for the XY 
model, where spin wave fluctuations do not renormal- 
ize the stiffness. In this limit, renormalization will only 
occur via vortex-antivortex fluctuations of the KT-flow 
equations, that eventually lead to Q. 



The solution to the flow equations (13) comes from 
noticing that f t [T(£) 2 / M A{£)] = 2T(£j z ^ 4 ' A{£), which 
allows us to implicitly construct the solution: 



T{£) 2 ' M A(£)e- 2i = C, 
(1 + A(£))e'^^ 2Ei -tt (mW)) 



A(£)T(£) 2 / M 



MT(£) 2 / M 




where C and D are constants determined from the bare 
parameters when £ = 0, and E n (x) is the exponential 
integral function. The RG flow is plotted in Fig. [T] for 
M = 2; other values of M look qualitatively the same. 



The RG flow equations ( 13 ) are for spin- waves only 



when the anisotropy A reaches a value of order 1 corre- 
sponding to the 0(2) phase, this must be supplemented 
by vortices and BKT arguments. We therefore stop the 
flow when A(^o) = 1 and ask the question: what is the 
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FIG. 2. Core size £o as a function of temperature T for vari- 
ous values of A. The red dashed line intersects the curves at 
Tbkt(A), showing the size of the vortex cores exactly at the 
unbinding transition. The lines end at the point where the 
vortex core size becomes comparable to the inter-vortex dis- 
tance, i.e. the point where £ no longer has any meaning and 
we are in the crossover from the 0(2) to the O(N) regime. 



renormalized value of temperature T(£q) at this scale? 
Above this scale, one sees only the physics of the tra- 
ditional XY model, and so T(£q) constitutes the initial 
condition of the BKT flow ([8| . As it is not even possible 
to define a vortex fugacity from the point of view of the 
0(N) spin waves, the previously discussed fugacity en- 
hancement no longer plays a role. In fact, from the point 
of view of the vortices in this approach, £o = ae l ° should 
be considered the microscopic (and only) length scale in 
the problem. Hence one finds a usual BKT transition at 
renormalized temperature T(£q) 
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to 



7r/2. Using Eqs 
trace this RG flow line back to the bare values of temper- 
ature and anisotropy, one finds that Tbkt ~ l/ m (l/A), 
in agreement with our result based entirely on vortices. 

We now discuss this length scale £ > which can of course 
be associated with the size of the vortex cores. Using 



( 14a I we have 



e 



a 
A 



T(e ) 

T 



2/M 



(15) 



Far below the BKT transition temperature, where T 
barely renormalizes, we can see that £o ~ a/V~A, again 
in accordance with our prior considerations. However, at 
higher temperatures T{£) flows toward strong coupling. 
Near Tbkt, where T(£q) ~ 1 as well, we obtain a loga- 
rithmic correction due to spin wave excitations that im- 
mediately leads to our result (JlJ. As the flow for £ < £q 
is governed by the RG equation of an 0(N) nonlinear 
(j-model, the vortex core enhancement is dictated by the 
high symmetry fixed point. In Fig. [2] we show the tem- 
perature dependence of the core size. 

Above but near Tbkt vortices proliferate, but the be- 
havior at large distances is still that of an 0(2) prob- 
lem. As temperature is further raised, a crossover occurs 
where A(£) never comes close to being of order unity, 
even as £ — > oo. A separatrix (shown in Fig. [T]) separates 



regions of flow where A(£ — > oo) does or does not di- 
verge, indicating the crossover from 0(2) to 0(N) behav- 
ior [22]. Analysing this behavior for small A yields that 
the crossover temperature also vanishes as l/ln(l/A), 
yet with a larger numerical coefficient: see Fig. [T] From 
the perspective of proliferated vortices, this crossover es- 
sentially corresponds to reaching the regime where the 
typical intervortex distance is of the order of the vortex 
core size, i.e. 0(N) fluctuations of the core govern the 
entire system. Following the BKT flow equation for the 
vortex fugacity confirms this interpretation. 

The full agreement between our two approaches is 
the most important conceptual conclusion of this work 
- above the length scale £oj the parameters flow with 
the usual BKT equations however below this length 
scale one can either choose to look at the problem from 



the point of view of vortices, (10), or from spin waves, 



( 13 ). In other words, it patches together the perturbative 



(spin- wave) and non-perturbative (vortex) aspects of the 
theory, something which has also recently been studied 
in a completely different context [2"5] . 

Finally we comment on a peculiarity of the case where 
M = 2, i.e. where we have two competing 0(2) order 
parameters. It has been established that this model at 
T = is governed by a tetracritical point, implying that 
both order parameters are non-zero at zero temperature 
for some range of A 24-26]. In this case we have two 
distinct BKT transitions for the two order parameters, 
where the upper transition takes place for the compo- 
nents that are stabilized by the anisotropy term A. A 
discussion of this special case will be given elsewhere [27] . 

In conclusion, we have analyzed a model of competing 
order parameters, at least one of which is of XY type. 
The energy balance between the competing states is con- 
trolled by a parameter A, such that the order param- 
eter symmetry is enhanced at A = 0. As the ground 
state of the system is fully ordered, we have a situa- 
tion where the ordering temperature vanishes as A — > 
without having quantum critical fluctuations: this van- 
ishing occurs solely because of the sensitivity with re- 
spect to spin waves of distinct order parameter sym- 
metries. More specifically, the Berezinskii-Kosterlitz- 
Thouless transition temperature Tbkt of vortex prolif- 
eration vanishes as A — > via a divergent core size, i.e. 
the integrity of topologically stable vortex configurations 
is undermined from within. Inside the giant core of such 
vortices, high symmetry spin-wave fluctuations further 
increase the core size. 

We believe that the observation of giant vortices and 
of intra-core excitations (see e.g. [35]) could be an im- 
portant clue in revealing the competing nature of order 
parameters in correlated many body systems. 
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